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1. Introduction 

In this paper we discuss eigenfunction expansions for a fundamental solution of Laplace's 
equation in the hyperboloid model of dimensional hyperbolic geometry. In particular, 
for a fixed R G (0, oo) and d >2, we derive and discuss Fourier cosine and Gegenbauer 
polynomial expansions in rotationally invariant coordinate systems, for a previously 
derived (see Cohl & Kalnins (2011) [8]) spherically symmetric fundamental solution 
of the Laplace-Beltrami operator in the hyperboloid model of hyperbolic geometry. 
Useful background material relevant for this paper can be found in Vilenkin (1968) |29j . 
Thurston (1997) [27], Lee (1997) [20] and Pogosyan & Winternitz (2002) [26]. 

This paper is organized as follows. In section [2], for the hyperboloid model of 
(i- dimensional hyperbolic geometry, we describe some of its global properties, such as 
its geodesic distance function, geodesic polar coordinates, Laplace-Beltrami operator 
(Laplacian), radial harmonics, and several previously derived equivalent expressions for 
a radial fundamental solution of Laplace's equation. In section ^ ioi d > 2, we derive 
and discuss Fourier cosine series for a fundamental solution of Laplace's equation about 
an appropriate azimuthal angle in rotationally invariant coordinate systems, and show 
how the resulting Fourier coefficients compare to the those in Euclidean space. In 
section m for d > 2, we compute Gegenbauer polynomial expansions in geodesic polar 
coordinates, for a fundamental solution of Laplace's equation in the hyperboloid model 
of hyperbolic geometry. In section [5] we discuss possible directions of research in this 
area. 

Throughout this paper we rely on the following definitions. For ai,a2,... G C, 
if z,j G Z and j < i then J2n=i^n = and nn=j '^'i ~ ^- "^^^ natural 
numbers is given by N := {1, 2, . . .}, the set Nq := {0, 1, 2, . . .} = N U {0}, and the set 
Z := {0, ±1, ±2, . . .}. The set R represents the real numbers. 

2. Global analysis on the hyperboloid 

2.1. The hyperboloid model of hyperbolic geometry 

Hyperbolic space, developed independently by Lobachevsky and Bolyai around 1830 
(see Trudeau (1987) [28]), is a fundamental example of a space exhibiting hyperbolic 
geometry. Hyperbolic geometry is analogous to Euclidean geometry, but such that 
Euclid's parallel postulate is no longer assumed to hold. There are several models of 
(i- dimensional hyperbolic space H^, including the Klein, Poincare, hyperboloid, upper- 
half space and hemisphere models (see Thurston (1997) [2^). The hyperboloid model 
for d-dimensional hyperbolic geometry (hereafter referred to as the hyperboloid model, 
or more simply, the hyperboloid), is closely related to the Klein and Poincare models: 
each can be obtained projectively from the others. The upper-half space and hemisphere 
models can be obtained from one another by inversions with the Poincare model (see 
section 2.2 in Thurston (1997) [27]). The model of hyperbolic geometry which we will 
be focusing on in this paper, is the hyperboloid model. 
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Minkowski space R*^'^ is a. {d+ l)-dimensional pseudo-Riemannian manifold which 
is a real finite-dimensional vector space, with Cartesian coordinates given by x = 
{xq,Xi, . . . ,Xd)- The hyperboloid model, also known as the Minkowski or Lorentz 
models, represents points in this space by the upper sheet (xq > 0) of a two- 
sheeted hyperboloid embedded in the Minkowski space R'^'^. It is equipped with a 
nondegenerate, symmetric bilinear form, the Minkowski bilinear form 

[x, y] = xoyo - Xiyi - ... - XdVd- 

The above bilinear form is symmetric, but not positive-definite, so it is not an inner 
product. It is defined analogously with the Euclidean inner product for R''"''^ 

(x, y) = Xoyo + + . . . + x^yd- 

The variety [x, x] = Xq — x\ — . . . — x^^ = R^, for x G R*^'^, using the language of 
Beltrami (1869) [3] (see also p. 504 in Vilenkin (1968) [29]), defines a pseudo-sphere of 
radius R. Points on the pseudo-sphere with zero radius coincide with a cone. Points on 
the pseudo-sphere with radius greater than zero lie within this cone, and points on the 
pseudo-sphere with purely imaginary radius lie outside the cone. 

For a fixed R G (0, oo), the i?-radius hyperboloid model is a maximally symmetric, 
simply connected, d-dimensional Riemannian manifold with negative- const ant sectional 
curvature (given by —l/R^, see for instance p. 148 in Lee (1997) [20]), whereas Euclidean 
space R*^ equipped with the Pythagorean norm, is a Riemannian manifold with zero 
sectional curvature. The hypersphere S'^, is an example of a space (submanifold) with 
positive-constant sectional curvature (given by 1/R^). 

In our discussion of a fundamental solution for Laplace's equation in the hj^erboloid 
model H^, we focus on the positive radius pseudo-sphere which can be parametrized 
through subgroup-type coordinates, i.e. those which correspond to a maximal subgroup 
chain 0{d,l) D ... (see for instance Pogosyan & Winternitz (2002) [26]). There exist 
separable coordinate systems which parametrize points on the positive radius pseudo- 
sphere (i.e. such as those which are analogous to parabolic coordinates, etc.) which can 
not be constructed using maximal subgroup chains (we will no longer discuss these). 

Geodesic polar coordinates are coordinates which correspond to the maximal 
subgroup chain given hj 0{d,l) D 0{d) D .... What we will refer to as standard 
geodesic polar coordinates correspond to the subgroup chain given by 0{d,l) D 0{d) D 
0{d — 1) D ■■■ D 0(2). Standard geodesic polar coordinates (see Olevskii (1950) 
[23] : Grosche, Pogosyan & Sissakian (1997) [IE]), similar to standard hyperspherical 
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coordinates in Euclidean space, can be given by 
Xq = R cosh r 
Xi = i?sinhrcos^i 
X2 = R sinh r sin 6i cos 62 



Xd-2 



(1) 



R sinh r sin 9i - ■ ■ cos 9d^2 

R sinh r sin ^1 ■ ■ ■ sin $^^2 cos 

i? sinh r sin ^1 ■ ■ ■ sin 6*^-2 sin 0, 

where re [0, cx)), G [0, 2it), and 6'i G [0, tt] for z G {1, . . . , ci - 2}. 

In order to study fundamental solutions on the hyperboloid, we need to describe 
how one computes distances in this space. One may naturally compare distances on the 
positive radius pseudo-sphere through analogy with the i?-radius hypersphere. Distances 
on the hypersphere are simply given by arc lengths, angles between two arbitrary 
vectors, from the origin, in the ambient Euclidean space. We consider the 0?- dimensional 
hypersphere embedded in R'^"'"^. Points on the hypersphere can be parametrized using 
hyperspherical coordinate systems. Any parametrization of the hypersphere S"^, must 



have 



x,x 



xi + 



+ 4 



R^, with R > 0. The distance between two points on the 



hypersphere is given by 

d(x, x') = i?7 



Rcos' 



X, X 



Rcos 



i?2 



(x, x'; 



(2) 



(x, x)(x',x' 

This is evident from the fact that the geodesies on S'^ are great circles (i.e. intersections 
of S'^ with planes through the origin) with constant speed parametrizations (see p. 82 
in Lee (1997) [20]). 

Accordingly, we now look at the geodesic distance function on the (i-dimensional 
positive radius pseudo-sphere H^. Distances between two points on the positive radius 
pseudo-sphere are given by the hyperangle between two arbitrary vectors, from the 
origin, in the ambient Minkowski space. Any parametrization of the hyperboloid H^, 
must have [x, x] = R^. The geodesic distance p G [0, 00) between two points x, x' G 
is given by 



(i(x, x') = i?cosh 



[x,x'] 



-Rcosh 



(3) 



[x,x][x',x'; 

where the inverse hyperbolic cosine with argument x G (1, 00) is given by (see (4.37.19) 
in Olver et al. (2010) [25]) cosh ~^ a: = log (x + — l) . Geodesies on are great 
hyperbolas (i.e. intersections of with planes through the origin) with constant speed 
parametrizations (see p. 84 in Lee (1997) [20]). We also define a global function 
p : H*^ X H'^ — [0, 00) which represents the projection of global geodesic distance 
function (|3]) on onto the corresponding unit radius hyperboloid H"^, namely 

p(x,X') :=rf(x,x')/i?, (4) 

where x = x/i? and x' = x'/i?. Note that when we refer to d{St, x') below, we specifically 
mean that projected distance given by dl]). 
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2.2. Laplace's equation and harmonics on the hyperboloid 

Parametrizations of a submanifold embedded in either a Euclidean or Minkowski space 
is given in terms of coordinate systems whose coordinates are curvihnear. These 
are coordinates based on some transformation that converts the standard Cartesian 
coordinates in the ambient space to a coordinate system with the same number of 
coordinates as the dimension of the submanifold in which the coordinate lines are curved. 

The Laplace-Beltrami operator (Laplacian) in curvilinear coordinates ^ = 
(^^,...,^'^) on a (i-dimensional Riemannian manifold (a manifold together with a 
Riemannian metric g) is given by 

where If?] = | det{gij)\, the infinitesimal distance is given by 

d 

d^' = E 9ijdCd^', (6) 



and 



i=l 

where Sj G {0, 1} with i,j G Z, is the Kronecker delta defined such that 
f 1 if « = ?, 

For a submanifold, the relation between the metric tensor in the ambient space and gij 
of dSD and (El) is 

, , >^ dx^ dx'' 

k,l=0 ^ ^ 

The ambient space for the hyperboloid is Minkowski, and therefore Gij = 
diag(l,-l,...,-l). 

The set of all geodesic polar coordinate systems corresponds to the many ways one 
can put coordinates on a hyperbolic hypersphere, i.e., the Riemannian submanifold 
U C defined for a fixed x' G such that (i(x, x') = 6 = const, where 
b G (0, oo). These are coordinate systems which correspond to subgroup chains starting 
with 0{d,l) D 0{d) D ■ ■ ■, with standard geodesic polar coordinates given by ([I]) being 
only one of them. (For a thorough description of these see section X.5 in Vilenkin (1968) 
[29].) They all share the property that they are described by ((i+l)-variables: r G [0, oo) 
plus d-angles each being given by the values [0,27r), [0,7r], [—11/2,11/2] or [0,7r/2] (see 
Izmest'ev et al. (1999, 2001) [HI EH]). 

In any of the geodesic polar coordinate systems, the geodesic distance between two 
points on the submanifold is given by (cf. ([3])) 

(i(x, x') = R cosh"'^ (cosh r cosh r' — sinh r sinh r' cos 7) , (8) 
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where 7 is the unique separation angle given in each hyperspherical coordinate system. 
For instance, the separation angle in standard geodesic polar coordinates ([1]) is given 
by the formula 

d-2 d-2 i-1 

COS 7 = cos(0 — (j)') JJ^ sin 6'jsin 9i + cos 6'jcos 6'/ JJ^ sin 6'jsin 6j' . (9) 

j=l i=l j=l 

Corresponding separation angle formulae for any geodesic polar coordinate system can 
be computed using ([2]), ([3]), and the associated formulae for the appropriate inner- 
products. 

The infinitesimal distance in a geodesic polar coordinate system on this submanifold 

is 

ds'^ = R^idr'^ + sinh^ r ^7^), (10) 

where an appropriate expression for 7 in a curvilinear coordinate system is given. If one 
combines ([I]), ([5]), ([9]) and (fTOl) . then in a particular geodesic polar coordinate system, 
Laplace's equation on is 

1 



d'^f df 1 

+ {d — 1) cothr— — I — — —As,d~if 



0, 11 



i9r2 dr sinh^ r 

where Agd-i is the corresponding Laplace-Beltrami operator on S'^"^ with unit radius. 

From this point onwards, S'^~^ will always refer to the {d — l)-dimensional unit 
hypersphere, which is a compact Riemannian submanifold with positive constant 
sectional curvature, embedded in R*^ and given by the variety + . . . + = 1. 

Geodesic polar coordinate systems partition into a family of concentric [d— 1)- 
dimensional hyperspheres, each with a radius r G (0, cxd), on which all possible 
hyperspherical coordinate systems for S*^"^ may be used (see for instance, in Vilenkin 
(1968) [29]). One then must also consider the limiting case for r = to fill out all of H^. 
In standard geodesic polar coordinates one can compute the normalized hyperspherical 
harmonics in this space by solving the Laplace equation using separation of variables 
which results in a general procedure which is given explicitly in Izmest'ev et al. (1999, 
2001) [171 [18]. These angular harmonics are given as general expressions involving 
trigonometric functions, Gegenbauer polynomials and Jacobi polynomials. 

The harmonics in geodesic polar coordinate systems are given in terms of a radial 
solution multiplied by the angular harmonics. The angular harmonics are eigenfunctions 
of the Laplace-Beltrami operator on S"'"^ with unit radius which satisfy the following 
eigenvalue problem 

As.-iy,^(x) = -/(/ + d- 2)F/^(x), (12) 

where x G S'^"^, Y;^(x) are normalized hyperspherical harmonics, / G Nq is the angular 
momentum quantum number, and K stands for the set of {d — 2)-quantum numbers 
identifying degenerate harmonics for each /. The degeneracy as a function of the 
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dimension d tells you how many linearly independent solutions exist for a particular 
/ value. The hyperspherical harmonics are normalized such that 



I' xK' 
K ) 



Y,^ {^)Yf {^)duj = 5\ 5\ 



where du is a volume measure on S''"^ which is invariant under the isometry group 
SO{d) (cf. (JH])), and for x + iy = z G C, ^ = x — iy, represents complex conjugation. 
The generalized Kronecker delta 6^' (cf. (I7j)) is defined such that it equals 1 if all of 
the {d — 2)-quantum numbers identifying degenerate harmonics for each I coincide, and 
equals zero otherwise. 

Since the angular solutions (hyperspherical harmonics) are well-known (see Chapter 
IX in Vilenkin (1968) [29]; Chapter 11 in Erdelyi et al. (1981) [I2]), we will now focus 
on the radial solutions, which satisfy the following ordinary differential equation 

d^u , , ^ , ^ du 1(1 + d — 2) 

— + {d-l) cothr- ^ ^-u = 0. 

ar"= dr sinh r 

Four solutions to this ordinary differential equation ufi_, ■ (loo) — t- C are given by 

%±(coshr) = . -Pd/2-r^^'^(cQshr), 
smh ' r ' 

and 



«2±(coshr) = . ,22-1 Qd/2-i ^^^^(coshr) 



where Pjf^Q'^ : (l,oo) — )■ C are associated Legendre functions of the first and second 
kind respectively (see for instance Chapter 14 in Olver et al. (2010) [25]). 

2.3. Fundamental solution of Laplace's equation on the hyperboloid 

Due to the fact that the space is homogeneous with respect to its isometry group, the 
pseudo-orthogonal group SO{d, 1), and therefore an isotropic manifold, we expect that 
there exist a fundamental solution on this space with spherically symmetric dependence. 
We specifically expect these solutions to be given in terms of associated Legendre 
functions of the second kind with argument given by cosh r. This associated Legendre 
function naturally fits our requirements because it is singular at r = and vanishes 
at infinity, whereas the associated Legendre functions of the first kind, with the same 
argument, are regular at r = and singular at infinity. 

In computing a fundamental solution of the Laplacian on H^, we know that 

-AH^(x,x') = 53(x,x'), (13) 

where g is the Riemannian metric on and 5g(x, x') is the Dirac delta function on 
the manifold H^. The Dirac delta function is defined for an open set U C with 
x, x' G such that 



(5g(x, x')dvo\g 



u 



if x' G f/, 
if x' ^ f/. 
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where dvolg is a volume measure, invariant under the isometry group SO{d, 1) of the 
Riemannian manifold H^, given (in standard geodesic polar coordinates) by 



dvoL = R'^ smh.'^-\duj :-- 



sinh^-^ r sin'^^^ 



■ sin 92d9i ■ ■ ■ d6, 



(14) 



Notice that as r — t- 0^ that dvolg goes to the Euclidean measure, invariant under 
the Euclidean motion group E{d), in spherical coordinates. Therefore in spherical 
coordinates, we have the following 



(5g(x,x'; 



9', 



R'^ sinh' 



sin ^2 ■ ■ ■ sin'* 



9r 



(15) 



d-l 



In general since we can add any harmonic function to a fundamental solution 
for the Laplacian and still have a fundamental solution, we will use this freedom 
to make our fundamental solution as simple as possible. It is reasonable to expect 
that there exists a particular spherically symmetric fundamental solution 'H^(x, x') on 
the hyperboloid with pure radial p(x, x') = d{'x,x.')/R (cf. (jlj)) and constant angular 
dependence (invariant under rotations centered about the origin), due to the influence of 
the point-like nature of the Dirac delta function. For a spherically symmetric solution to 
the Laplace equation, the corresponding Agd-i term vanishes since only the / = term 
survives. In other words, we expect there to exist a fundamental solution of Laplace's 
equation such that ^{^(x, x') = /(p). 

In Cohl & Kalnins (2011) [8], we have proven that on the -R-radius hyperboloid 
H^, a fundamental solution of Laplace's equation can be given as follows. 

Theorem 2.1. Let de {2,3,.. .}. Define Id : (0, oo) ^ R as 

/"°° dr 

Up)--- 



sinh*^ ^ X ' 



X, x' G H^, and Uj^ 



[H^ X H^) \ {(x, x) : X e H^} -> R defined such that 



T-Ljii^, X 



r (rf/2) 



where p := cosh"^ ([x, x']) is the geodesic distance between^ andx' on the pseudo-sphere 
of unit radius H^, with x = x/i?, x' = x'/i?, then "H^ is a fundamental solution for — A 
where A is the Laplace- Beltrami operator on H^. Moreover, 



Up) 



{-if 



{d- 

(_l)('i-l)/2 



2)!! 
'id 



log coth ^ + cosh p 



d/2-l 



k=l 



(2A;-2)!!(-l)'= 
(2A:- l)!!sinh2^p 



if d even, 



3)!! 



_(rf-2)!! 

3\, ^ 



-1)'' coth 



2k-l 



k=l 



{2k-l){k-l)\{{d-2k-l)/2)\ 



} if d odd. 



or 



(_l)('^-l)/2 



(d-3)! 
Jd^l 



(d-l)/2 

1 + cosh p 



k=l 



(2fc-3)!!(-l)fe 
(2A;-2)!!sinh2'=-\ 
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(rf - 1) cosh'^^ V V 2 '2' 2 cosh> 



((i — 1) coshp sinh'^ \2' ' 2 'cosh^p/' 

^ 2'^/2-ir (rf/2) sinh-^/^-i p Q./2-i(coshp), 

where !! ^s i/ie double factorial, 2F1 is the Gauss hypergeometric function, and Q'^ is the 
associated Legendre function of the second kind. 

For a proof of this theorem, see Cohl & Kalnins (2011) [S]. 

3. Fourier expansions for a Green's function on the hyperboloid 

Now we compute the Fourier expansions for a fundamental solution of the Laplace- 
Beltrami operator on H^. 

3.1. Fourier expansion for a fundamental solution of the Laplacian on H|, 

The generating function for Chebyshev polynomials of the first kind (Fox & Parker 
(1968) [11], p. 51) is given as 

2 °° 

i + (16) 

n=0 

where \z\ < 1, r„ : [—1, 1] — R is the Chebyshev polynomial of the first kind defined as 
Ti{x) := cos(/ cos""^ x), and e„ := 2 — 5^ is the Neumann factor (see p. 744 in Morse & 
Feshbach (1953) [22]), commonly-occurring in Fourier cosine series. If substitute z = 
with rj G (0, 00) in ( lT6l) . then we obtain 

sinh?7 



cosh rj — cos 



^e„cos(n^/;)e-"''. (17) 

n=0 

Integrating both sides of f|T7j) with respect to 77, we obtain the following formula 
(cf. Magnus, Oberhettinger & Soni (1966) [21], p. 259) 



log {1 + z^ - 2z cos tfj) = -2Y^ ^—^z''. 



n 

n=l 



In Euclidean space R"', a Green's function for Laplace's equation (fundamental 
solution for the Laplacian) is well-known and is given in the following theorem (see 
Folland (1976) [I3]; p. 94, Gilbarg & Trudinger (1983) [I5]; p. 17, Bers et al. (1964) [1], 
p. 211). 
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Theorem 3.1. Let d eN. Define 



rW2) „^_^/||2-d ifrf=lorrf>3, 



— log||x-x'||"^ ifd = 2, 

27r 

then Q'^ is a fundamental solution for — A in Euclidean space 'Rf' , where A is the Laplace 
operator in R°'. 

Therefore if we take z = r^/r^ in ( |T8l) . where := ^^^{'^,1^'} with r, r' G [0,cx3), 
then using polar coordinates, we can derive the Fourier expansion for a fundamental 
solution of the Laplacian in Euclidean space for d = 2 (cf. Theorem 13. ip . namely 

2 1 II /|| 1 -A cos(n(0- 0')) /r<\'^ 

:= log x-x = logr> - > — , (19) 

^-^ n Vr>/ 

n=l ^ 

where = —27iQ'^ (cf. Theorem 13. ip . On the hyperboloid for d = 2 we have a 
fundamental solution of Laplace's equation given by 

r 2 T .1 1 1 1 cosh d{% X') + 1 

" ■= ""^ S"""' ^ ) = 5 coshd(x.x-)-l ' 

where 1)^ = 27rn% (cf. Theorem O and (EI]) below). Note that because of the R"^-^ 
dependence of a fundamental solution of Laplace's equation for c? = 2 in Theorem 12. 1[ 
there is no strict dependence on R for T-L^ or f)^, but will retain the notation nonetheless. 
In standard geodesic polar coordinates on H|. (cf. ([I])), using and cos 7 = cos(0 — 0') 
(cf. ([9])) produces 

cosh (i(x, x') = cosh r cosh r' — sinh r sinh r' cos(0 — </>'), 

therefore 

9 1, cosh r cosh r' + 1 — sinh r sinh r' cos((7!) — 0') 

f) = - log — —. 

2 cosh r cosh r' — 1 — sinh r sinh r' cos(<^ — (j)') 

Replacing tp = (p — (p' and rearranging the logarithms yield 

,2 1, cosh r cosh r' + 1 1 , ,n 1, /. ,n 

^ = 2 coshrcoshr--l + 2 " ^+ ^) - - log (1 - cos ^) , 

where 

sinh r sinh r' 
^ cosh r cosh r' ± 1 
Note that 2± G (0, 1) for r, r' G (0, 00). We have the following MacLaurin series 

°° 

log(l -x) = - V— , 

n 

n=l 



where a: G [— 1, 1). Therefore away from the singularity at x = x' we have 

00 ^ 

A± := log(l - 2±cosV^) = -^^cos^V- (20) 



k 

k=l 
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We can expand the powers of cosine using the following trigonometric identity 

k 

i 1 

COS" 



2k 

n=0 

which is the standard expansion for powers using Chebyshev polynomials (see for 
instance p. 52 in Fox & Parker (1968) [S]). Inserting this expression in fl^Uj) . we obtain 
the following double-summation expression 

oo k k f 1 \ 

A± = -5^5^^| lcos[(2n-fc)^]. (21) 

k=l n=0 



Now we perform a double- index replacement in ( 12T]) . We break this sum into two 
separate sums, one for k < 2n and another for k > 2n. There is an overlap when both 
sums satisfy the equality, and in that situation we must halve after we sum over both 
sums. If /c < 2n, make the substitution k' = k — n and n' = 2n — k. It follows that 
k = 2k' + n' and n = n' + k', therefore 

k\ _ f2k' + n'\ _ hk' + n' 
nj~yn' + k'j~yn' + k' 

\i k>2n make the substitution k' = n and n' = k — 2n. Then k = 2k' + n' and n = k', 
therefore 

'k\ _ f 2k' + n'\ _ f 2k' + n' 
n I \ n j \ k' + n' 

where the equalities of the binomial coefficients are confirmed using the following identity 



n 



n 



k 



where n, G Z, except where /c < or n — /c < 0. To take into account the 
counting which occurs when k = 2n (which occurs when n' = 0), we introduce 
of e„'/2 into the expression (and relabel k' ^ k and n' ^ n). We are left with 




A. 



If we substitute 



2 cos{nip) 



^2k+n 



n=l 



k=0 




22^-+" (2 A; 



double- 
a factor 



(22) 



k\ 



into the first term of (l22l) . then we obtain 




^ oo 
k=l 



2k 



k\k 

We are left with 
/± = -log2-|-log I 1 



dz'^ ^ (2)^. 



, 2k 



k=l 



k\ 



"± dz'^ 



l-z: 



1 2 



log2 + log / (coshr>±l)(coshr< + l) ' 
\ cosh r cosh r' ± 1 
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12 



-,2fc 




n 



1 



n 



into the second term of (1221), then the Fourier coefficient reduces to 



J. 



± 



>n-l /L^ 



n+1 



n 



2k+n 



k=0 



k\(n + l)k 



2k + n 



>n-l 



z± 



n+1 



n 



k=0 



k\{n 



kj 2/= 



The series in the integrand is a Gauss hypergeometric function which can be given as 



oo 

E 

A:=0 



n + 1 



n + 2 



k\(n 



k^2k 



2'^r?' / ^ 



where Pg " is an associated Legendre function of the first kind with vanishing degree 
and order given by —n. This is a consequence of 



2F1 ( a,b;a + b 



X 



-3/2p 



p:itxi2 (v^) 



2j v^r^ 

where x G (0, 1) (see for instance Magnus, Oberhettinger & Soni (1966) [21], p. 53), and 
the Legendre function is evaluated using (cf. (8.1.2) in Abramowitz & Stegun (1972) 



P,-^{x) 



1 



X 



X 



n/2 



where n G Nq. Therefore the Fourier coefficient is given by 



J± 

Finally we have 
A± 



yTZ4 1-42 ^1 + ^/ 



1 



n/2 



2 

n 



1 



^2 



-1 n/2 



log 2 + log 



/ (coshr> ± l)(coshr<; + 1) 



-2E 



cosh r cosh r' ± 1 
(coshr> =p l)(coshr< 



n=l 



n 



n/2 



(coshr> ± l)(coshr< + 1) J ' 

and the Fourier expansion for a fundamental solution of Laplace's equation for the d = 2 
hyperboloid is given by 



cosh r> + 1 
coshr 



> 



^ ^ cos(n(0-0O) 

n=l 



n 



coshr< 



coshr< + 1 



n/2 



coshr> + 1 
coshr> — 1 



n/2 



coshr- 



cosh 



n/2" 



■ (23) 



This exactly matches up to the Euclidean Fourier expansion (ITU]) as r, r' — )■ 0"*". 
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3.2. Fourier expansion for a fundamental solution of the Laplacian on H|j 

The Fourier expansion for a fundamental solution of the Laplacian in three- 
dimensional Euclidean space (here given in standard spherical coordinates x = 
(rsin^cos0, rsin^sin0, r cos^)) is given by (cf. Theorem 13.11 and see (1.3) in Cohl 
et al. (2001) [9]) 
1 



^ y e-(*-^')g^ ( r' + r"-2rr'cos9cos9'\ 

TrVrr'sin^sin^' ^ 1/2 2rr' sin ^ sin J 



m=— 00 



These associated Legendre functions, toroidal harmonics, are given in terms of complete 
elliptic integrals of the first and second kind (cf. (22-26) in Cohl & Tohline (1999) [TO]). 
Since (5-1/2(2;) is given through (cf. (8.13.3) in Abramowitz & Stegun (1972) [1]) 



the m = component for q is given by 



o, 2 ( I 4rr' sin 6^ sin 6^' , , 




^^^2 + ^f'i _ 2rr' cos{e + 6') \]l + r'^ - 2rr' cos{ 

A fundamental solution of the Laplacian in standard geodesic polar coordinates on 
Hfj is given by (cf. Theorem 12.11 and (13T1) below). 

cosh (i(x, x') 



:= cothci(x,x') - 1 



cosh^ (i(x, x') — 1 
cosh r cosh r' — sinh r sinh r' cos 7 



•y/ (cosh r cosh r' — sinh r sinh r' cos 7)^ — 1 

where f)^ = AirRH^, and x, x' G H^, such that x = x/i? and x' = x'/i?. In standard 
geodesic polar coordinates (cf. ([9])) we have 

cos 7 = cos 6* cos 6*' + sin ^ sin 6*' cos (0 — </)'). (25) 

Replacing ip = cj) — cj)' and defining 

A : = cosh r cosh r' — sinh r sinh r' cos ^ cos 6', 

and 

S := sinh r sinh r' sin ^ sin 6', 

we have in the standard manner, the Fourier coefficients HH"^ : [0, 00)^ x [0, tt]^ — ?> R of 
the expansion (cf. (132|1 below) 

00 

f)3(x,x') = 5^cos(m(0-0'))H^/'(r,r',^,^'), (26) 

m=0 
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defined by 



TT 



(cos^-^) (cosV^-^) 
If we make tlie substitution x = cosz/^, tliis integral can be converted to 
Hll^r, r', 6, 9') = -5° + ^ (A/E - x) T^(x)rfx _ 



TX 



^(l-x)(l+a:) (x-^) [x-^) 



wliere is tlie Cliebysliev polynomial of the first kind. Since Tm{x) is expressible as 
a finite sum over powers of x, (!28|) involves the square root of a quartic multiplied by a 
rational function of x, which by definition is an elliptic integral (see for instance Byrd 
& Friedman (1954) [5]). We can directly compute (128|) using Byrd & Friedman (1954) 
(0, (253.11)). If we define 

A-l A+l 
d := -1, y := -1, c := 1, 6 := — a := — (29) 

(clearly d<y<c<b<a), then we can express the Fourier coefficient (l28ll . as a linear 
combination of integrals, each of the form (see Byrd & Friedman (1954) [5], (253.11)) 



x^dx 



d'g 





'1 


2 2 1 

— a^sn u 


Jo 


1 


— a^sn^M 



p 



du, (30) 



A/(a — x)(6 — x)(c — x)(x — d) 

where p G {0, . . . ,m + 1}. In this expression sn is a Jacobi elliptic function (see for 
instance Chapter 22 in Olver et al. (2010) [25]). 

Byrd & Friedman (1954) |5] give a procedure for computing ( !30l) for all m G Nq. 
These integrals will be given in terms of complete elliptic integrals of the first three 
kinds (see the discussion in Byrd & Friedman (1954) [5], p. 201, 204, and p. 205). To 
this effect, we have the following definitions from (253.11) in Byrd & Friedman (1954) 
[5], namely 

n C — d 

2 ^ b{c - d) 
c{b-d)' 
2 



^{a-c){b-dy 



^b-d)ic-y) 
^ic-d)ib-yy 

Ul = F{(f,k), 

2 ^ {a-b){c-d) 
{a-c){b-dy 

with k"^ < a^. For our specific choices in (129|) . these reduce to 
2 _ 25 
" ~ ^ + 5-1' 
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,2_ 2(A-1) 



a 



^ A+B-V 

2B 



^{A + B-1){A-B + 1) 
AB 

{A + B-1){A-B + 1)' 



and 



Ml = K{k). 
Specific cases include 



gKik) 



\J {a — x){b — x){c — x){x — d) 
(Byrd & Friedman (1954) [5], (340.00)) and 

xdx eg 



[alK{k) + (a^ - alma, k)] 

^/{a - x){b - x){c- x){x - d) «^ 



(Byrd & Friedman (1954) [5], (340.01)). 
In general we have 



x^dx c^ga^ipl ^ (a^ — a\y 

~ 2^ 7:27-77-— TT^i 



^{a-x){h-x){c-x){x-d) ^ ^2,^.,^^ _ 

(Byrd & Friedman (1954) [5], (340.04)), where 
K) = i^(A;), 
Vi = n(a,fc), 

^2 = 2 ^ [(^' - «')^(^) + «'^(^) + (2«'^' + 2«' - - 3A:2)n(a, k)] 

and larger values of Vj can be computed using the following recurrence relation 

1 

^""^^ " 2(m + 2)(l-a2)(A;2_a2) 

X [(2m + l)fcV„, + 2(m + l)(a^fc^ + a^ - ?>k^)V„,+i 

+ (2m + 3)(a^ - 2a^k'^ - 2a^ + ?>k'^)Vm+2] 
(see Byrd & Friedman (1954) [5], (336.00-03)). For instance, 

x'^dx 



y^(a:^^"xy(6^^~ry(cF^^"xy(:r^^^ 
2 

= ^ \a\K{k) + 2a2(c,2 _ a2)n(a, k) + (a^ _ ^2)2^1 
a^ 



(see Byrd & Friedman (1954) 0, (340.02)). 
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In general, the Fourier coefficients for i)^ will be given in terms of complete elliptic 
integrals of the ffist three kinds. Let's directly compute the m = component, in which 
( l28l) reduces to 

{A/B - x) dx 



(r, r', 9,9') = -! + -[' 



;)(l + x) (x-^) (x-^) 
Therefore using the above formulae, we have 
Um=o = \^T{ry,e,9') 

2K{k) 2{A - B -l)U{a,k) 



7r 



y/{A-B + l){A + B-l) 7iy/{A-B + l){A + B -1) 



2 

= — 1 H — {K{k) + [cosh r cosh r' — sinhr sinhr'cos(^ — 9') — 1] Il{a, k)} 

TT 

X [cosh r cosh r' — sinh r sinh r' cos{9 — 9') + 1] 

X [cosh r cosh r' — sinh r sinh r' cos{9 + 9') — 1] . 
Note that the Fourier coefficients 

I) |m=0 ~^ \m=0, 

in the limit as r, r' — 0+, where g^|m=o is given in ( l24l) . This is expected since is a 
manifold. 

3.3. Fourier expansion for a fundamental solution of the Laplacian on 

For the (i- dimensional Riemannian manifold H^, with d > 2, one can expand a 
fundamental solution of the Laplace-Beltrami operator in an azimuthal Fourier series. 
One may Fourier expand, in terms of the azimuthal coordinate, a fundamental solution 
of the Laplace-Beltrami operator in any rotationally-invariant coordinate systems 
which admits solutions via separation of variables. In Euclidean space, there exist 
non-subgroup-type rotationally invariant coordinate systems which are separable for 
Laplace's equation. All separable coordinate systems for Laplace's equation in d- 
dimensional Euclidean space R'^ are known. In fact, this is also true for separable 
coordinate systems on (see Kalnins (1986) [I9]). There has been considerable work 
done in two and three dimensions, however there still remains a lot of work to be done 
for a detailed analysis of fundamental solutions. 

We define an unnormalized fundamental solution of Laplace's equation on the unit 
hyperboloid t)'^ : (H'^ x H'^) \ {(x, x) : x G H"^} R such that 

[)^(x,x') := J,(p(x,x')) = -YJJJ^'^iM- (31) 

In our current azimuthal Fourier analysis, we therefore will focus on the relatively easier 
case of separable subgroup-type coordinate systems on H^, and specifically for geodesic 
polar coordinates. In these coordinates the Riemannian metric is given by (fTOj) and we 
further restrict our attention by adopting standard geodesic polar coordinates ([T]). 
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In these coordinates would would like to expand a fundamental solution of Laplace's 
equation on the hyperboloid in an azimuthal Fourier series, namely 



f)^(x, x') = J2 cos(m(0 - (j)')^'-' (r, r', ^i, . . . , ^,_2, e[, . . . , 6',^,) (32) 

m=0 

where Hm^^^ : [0, oo)^ x [0, vr]^"'"'^ — R is defined such that 

^T-\r, r',e,,..., e,.2, e[,..., := — r x') cos(m(</) - <j>'))d{<jy - 0') (33) 

Jo 

(see for instance Cohl & Tohline (1999) [10]). According to Theorem 12.11 and fl3Tl) . we 
may write f)'^(x, x') in terms of associated Legendre functions of the second kind as 
follows 



- ) = 2^/^-ir(d/2)(sinhrf(x,xO)^/^-^ ^^/^^^ ^"""^ "'^^ ■ ^^'^ 
By ([3]) we know that in any geodesic polar coordinate system 

cosh (i(x, x') = cosh r cosh r' — sinh r sinh r' cos 7, (35) 

and therefore through fl33l) . fl34l) . and fl35l) . in standard geodesic polar coordinates, the 
azimuthal Fourier coefficient can be given by 

H'J,'-\r,r\e,,...,e,.,,e[,...,e',_,) 

-in{d/2~i) Q'^Jjl~\ {A-B COS ^) cos(mV;) (36) 



-iTT(d/2-l) 



-dip, 



2'^/2-ivrr(rf/2) J, [(A - 5COS V')^ - 1]^'-'^/' 
where ijj := (j) - (p' , A, B : [0, oo)^ x [0, vr]^"'-^ ^ R are defined through and ((351) as 

y4(r, r',9i, . . . , 9^-2, 9[, . . . , 0'^_2) := cosh r cosh r' cos OiCosO/ Y\ sin 6'jsin6'/, 

i=i j=i 

and 

d-2 

B{r, r',9i, . . . , 6^-2, d'l, ■ ■ ■ , d'd~2) '■— sinh r sinh r' JJ^ sin6'jsin^/. 

i=l 

Even though (l36l) is a compact expression for the Fourier coefficient of a 
fundamental solution of Laplace's equation on for d G {2,3,4,...}, it may be 
informative to use any of the representations of a fundamental solution of the Laplacian 
on from Theorem 12. II to express the Fourier coefficients. For instance if one uses the 
finite-summation expression in the odd-dimensions, on can write the Fourier coefficients 
as a linear combination of integrals of the form 

[{a + h)/2-xf''~\Pdx 

1 (a — xY~^{h — xY^^ \J {a — x){b — x) (c — x){x — d) 

where x = cosip, k G {l,...,{d— l)/2}, p G {0, ...,m}, and we have used the 
nomenclature of section 13. 2[ This integral is a rational function of x multiplied by 
an inverse square-root of a quartic in x. Because of this and due to the limits of 
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integration, we see that by definition, these are all given in terms of complete elliptic 
integrals. The special functions which represent the azimuthal Fourier coefficients on 
are unlike the odd-half-integer degree, integer-order, associated Legendre functions 
of the second kind which appear in Euclidean space R*^ for d odd (see Cohl (2010) [6]; 
Cohl & Dominici (2010) [7j), in that they include complete elliptic integrals of the third 
kind (in addition to complete elliptic integrals of the first and second kind) (cf. section 
13. 2p in their basis functions. For d > 2, through (4.1) in Cohl & Dominici (2010) [7] 
and that is a manifold (and therefore must locally represent Euclidean space), the 
functions Hm^~^ are generalizations of associated Legendre functions of the second kind 
with odd-half-integer degree and order given by either an odd-half-integer or an integer. 



4. Gegenbauer expansion in geodesic polar coordinates 

In this section we derive an eigenfunction expansion for a fundamental solution of 
Laplace's equation on the hyperboloid in geodesic polar coordinates for c? G {3,4, . . .}. 
Since the spherical harmonics for d = 2 are just trigonometric functions with argument 
given in terms of the azimuthal angle, this case has already been covered in section [3?T1 

In geodesic polar coordinates, Laplace's equation is given by (cf. (fTTj) ) 



d^f df 1 



/ = 0, (37) 



^ dr"^ dr sinh r 

where / : R and Agd-i is the corresponding Laplace-Beltrami operator on the 

(d — l)-dimensional unit sphere S'^"^. Eigenfunctions : S'^"^ — )■ C of the Laplace- 
Beltrami operator Agd-i, where / G Nq and K is a set of quantum numbers which 
label representations for / in separable subgroup type coordinate systems on S'^"^ (i.e. 
angular momentum type quantum numbers, see Izmest'ev et al. (2001) [H]), are given 
by solutions to the eigenvalue problem ( TT^ . 

In standard geodesic polar coordinates ([1]), K = {ki, . . . , kas, \kd-2\) G Nq~^ with 
ko = I > ki > . . . > kd^3 > \kd-2\ > 0, and in particular kd-2 e {-kd-3, • • • , kd-a}- A 
positive fundamental solution T-Lj^ : (H^ x H^) \ {(x, x) : x G H^} — )■ R on the -R-radius 
hyperboloid satisfies ( IT3|) . The completeness relation for hyperspherical harmonics in 
standard hyperspherical coordinates is given by 

\"\"v^(f) R \vKm> — - ^^^^ ~ ^'^^ • • • ^^^'^-^ ~ ^^"1^ 

l^l^^i [Ou...,ed-i)Yi . . . , J - — -d-^n, — ' 

where K = [ki, . . . , kd-2) and I = k^ E Nq. Therefore through f|T5|) . we can write 



00 



For fixed r, r' G [0, 00) and 0^, . . . , 9'^_^ G [0, tt], since "H^ is harmonic on its domain, 
its restriction is in C^(S'^~^), and therefore has a unique expansion in hyperspherical 
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harmonics, namely 

oo 

nU^, ^') = Y.Y. ^^(^' ^i' • • • ' • • • ' ^-^-i)' (39) 

1=0 K 

where uf : [0, oof x [0, Tif^^ C. If we substitute dSHD and ([39]) into (fT3]) and use 
(I37j) and dn]), we obtain 

^ /(/ + rf_2)" 



, {d- l)cothr- . 
c(r^ dr smh r 



This indicates that for ui : [0, oo)^ — )■ R, 

«r (r, r', . . . , = r')F,^K,...,^:,_i), (41) 

and from (!39|) the expression for a fundamental of the Laplace-Beltrami operator in 
hyperspherical coordinates on the hyperboloid is given by 

oo 

<(x, x') = ui{r, r') J2 Yi^iO,, • • • , 9 ,^,)Y,^ {e[, . . . , 6',^,) . (42) 

1=0 K 

The above expression can be simplified using the addition theorem for hyperspherical 
harmonics (see for instance Wen & Avery (1985) [30], section 10.2.1 in Fano & Rau 
(1996), Chapter 9 in Andrews, Askey & Roy (1999) [2] and especially Chapter XI in 
Erdelyi et al. Vol. II (1981) [I2]), which is given by 

^ y,^(x) Wo = ^ Ji^fi (2/ + - 2)Cf^-^(cos7), (43) 



K 



where 7 is the angle between two arbitrary vectors x, x' G S'^^^ given in terms of ([2]). 
The Gegenbauer polynomials : [—1, 1] — R, / G Nq, Re/i > —1/2, can be defined in 
terms of the Gauss hypergeometric function as 



Cr(^) :=^2Fi(-/,/ + 2/i;/i 
The above expression ([42]) can be simplified using ([43]) . therefore 

^«(^'^') = 2Ji(f-9^ f;nKr,rO(2/ + rf-2)Cf '^(cos7). (44) 



=0 



Now we compute the exact expression for ui{r,r'). By separating the angular 
dependence in MOt) and using ([^T]) . we obtain the differential equation 

f^ + (d-l) cothr^ - Kl + d-2)ui ^ 6{r - r') 

dr'^ ^ ' dr sinhV i?'^-2 sinh'^" V' ' ^ ' 

Away from r = r', solutions to the differential equation f [lS]) must be given by solutions 
to the homogeneous equation, which are are given in section 12.21 Therefore, the solution 
to ( 1^5]) is given by 



n,(r, r') = P;ff-^+') (cosh r<)g;;^:;+' (cosh r>), (46) 

smh r smh r' 
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such that Ui{r,r') is continuous at r = r', where A G R. 

In order to determine the constant A, we first make the substitution 

Vi{r,r') = (sinhr sinhr')*''^~^-'''^u;(r, r'). (47) 

This converts (145|) into the following differential equation 



d'^Vi{r,r') 1 



(d - l + 2/)(d-3 + 2/) 



{d-1 



i2 



Mr,r) = j^j^. 



sinh^ r 

which we then integrate over r from r' — e to r' + e, and take the limit as e — 0+. We 
are left with a discontinuity condition for the derivative of vi (r, r') with respect to r 
evaluated at r = r', namely 

dvi{r, r') 



lim 

r-)-o+ dr 



(4J 



After inserting ( l46l) with (H71) into (l48i) . substituting 2; = coshr', evaluating at r = r', 
and making use of the Wronskian relation (e.g. p. 165 in Magnus, Oberhettinger & Soni 
(1966) [21]) 

which is equivalent to 



W { P-'^ (cosh r'),Q'^ (cosh r') } 



sinh r' 



we obtain 



and hence 



g-j7r(d/2-l+0 



g-i7r(d/2-l+0 

"'^^'^'^ = ^^^^(^hn^T;^^ 

and therefore through fl44p . we have 

, ^ r(ci/2) e-W^-^) 

' ^ 27r'^/2i?rf"2(d-2) (sinhrsinhr')'^/'-! 

00 

X 5^(-l)'(2/ + - 2)P,7(l/J"^+')(coshr<)Q2/^:;+'(coshr>)C'f/^-^(cos7). (49) 
/=o 

As an alternative check of our derivation, we can do the asymptotics for the 

product of associated Legendre functions P^2'^^^~^'''^'*(coshr<)(5^/2-i''''('^°^'^'^>) POl) 
as r, r' — )■ 0"*". The appropriate asymptotic expressions for P and Q respectively can be 
found on p. 171 and p. 173 in Olver (1997) [21]. For the associated Legendre function 
of the first kind there is 

" * ' r(^ + i) ' 
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as 2; — )■ 1, /i 7^ —1, —2, . . ., and for the associated Legendre function of the second kind 
there is 

^"^^ 2[(.-l)/2f/^' 

as z — )■ 1"*", Re /i > 0, and z/+/i 7^ —1, —2, —3, .... To second order the hyperbohc cosine 
is given by coshr ~ 1 +r^/2. Therefore to lowest order we can insert coshr< ~ 1 +r^/2 
and coshr> ~ 1 + r|/2 into the above expressions yielding 

(r</2)'^/2-i+' 



and 



2(r>/2)^/2-i+' 

as r,r' — )■ 0^. Therefore the asymptotics for the product of associated Legendre 
functions in (1491) is given by 

„i7r(d/2-l+0 \ l.+d/2-l 

r,)Q*-r'(cosh r,) ~ i^-^ (li) (50) 

(the factor 2/ + (i — 2 is a term which one encounters regularly with hyperspherical 
harmonics). Gegenbauer polynomials obey the following generating function 

^ oo 

J2Ct{x)z\ (51) 



(1 + ^2 _ 2zxf 

where x G [—1, 1] and \z\ < 1 (see for instance, p. 222 in Magnus, Oberhettinger & Soni 
(1966) [21]). The generating function for Gegenbauer polynomials (15T!) can be used to 
expand a fundamental solution of Laplace's equation in Euclidean space R"^ (for d > 3, 
cf. Theorem 13. ip in hyperspherical coordinates, namely 



oo 



J 



||,_,.||.-. = E 3^<'"'(-^7), (52) 

" " 1=0 > 

where 7 was defined in fH3|) . Using flS2]) and Theorem 13. ![ since — )■ 
sinhr, sinhr' — )■ r, r' and (150!) is satisfied to lowest order as r, r' — )■ O"*", we see that (149!) 
obeys the correct asymptotics and our fundamental solution expansion locally reduces 
to the appropriate expansion for Euclidean space, as it should since is a manifold. 

Note that ( H9|l can be further expanded over the remaining (ci— 2)-quantum numbers 
in K in terms of a simply separable product of normalized harmonics Y';^(x)y5^(x'), 
where x, x' G S'^"^, using the addition theorem for hyperspherical harmonics ( 143|) (see 
Cohl (2010) [0] for several examples). 

It is intriguing to observe how one might obtain the Fourier expansion for d = 2 
(125]) from the the expansion (jinD, which is strictly valid for d > 3. If one makes 
the substitution fi = d/2 — 1 in (149|) then we obtain the following conjecture (which 
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matches up to the generating function for Gegenbauer polynomials in the Euclidean 
hmit r, r' — )■ 0'*') 



1 



sinh'* p 



Q'icoshp) 



2^r(/x + 1) 

(sinhr sinhr')/^ 



X 



^(_l)n!i±i^p^;(M+n)(^^gl^^^)gM+n(coshr>)C„^(cOS7), (53) 



n=0 



for all /i G C such that Re^u > — 1/2. If we take the limit as /i — )■ in fl55]) and use 

enTnix) (54) 



liin !?-L^C/^(a;) 



(see for instance (6.4.13) in Andrews, Askey & Roy (1999) [2]), where T„ : [—1, 1] — >■ R 
is the Chebyshev polynomial of the first kind defined as T;(x) := cos(/ cos"^ x), then we 
obtain the following formula 



cosh p + 1 



- log """"^ — = 

2 cosh p — 1 ^-^ 

n=0 



cosh p — 

where cosh p = cosh r cosh r' — sinh r sinh r' cos 
following formulae 



l)"Po~"(coshr<)g;j(coshr>) cos(n(0 - 0')), 

- 0'). By taking advantage of the 



n-"(^) = ^ 
n! 



for n > 0, 



z - 1 

z+ 1, 

Qo(.) = 2log73T 



n/2 



(55) 



(56) 



(57) 



for n > 1 then fl25]) is reproduced. The representation fl55l) follows easily from the Gauss 
hypergeometric representation of the associated Legendre function of the first kind (see 
(8.1.2) in Abramowitz & Stegun (1972) P) 



(see (8.4.2) in Abramowitz & Stegun (1972) OP), and 



Q^(^) = -(-l)"(n-l)! 





'z+l 


n/2 


'z - 1" 






_z-l_ 




.^ + 1. 





P,f(z) 



z + l 



r(i -/i) [z - 1 



m/2 



2-ri 



-z/, 1/ + 1; 1 - p; 



1-z 



{5t 



One way to derive the representation of the associated Legendre function of second 
kind (!57|) is to use Whipple formula for associated Legendre functions (cf. (8.2.7) in 
Abramowitz & Stegun (1972) [1]) 



0F2r(z/ + p + l){z^- l)-l/4g»-/.p— 1/2 



and (8.6.9) in Abramowitz & Stegun (1972) [T], namely 



2f2^ 



1) 



-1/4 



TT 2Z/+ 1 



+ V2 



1 



n i^+1/2 



+ 72 



-1/-1/2 



for ^ -1/2. 
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4.I. Addition theorem for the azimuthal Fourier coefficient on H|j 

One can compute addition theorems for the azimuthal Fourier coefficients of a 
fundamental solution for Laplace's equation on for > 3 by relating directly 
obtained Fourier coefficients to the expansion over hyperspherical harmonics for the 
same fundamental solution. By using the expansion of 'H^(x, x') in terms of Gegenbauer 
polynomials fH9|) in combination with the addition theorem for hyperspherical harmonics 
( l43l) expressed in, for instance, one of Vilenkin's polyspherical coordinates (see section 
IX.5.2 in Vilenkin (1968) [29]; Izmest'ev et al. (1999,2001) [HI [18]), one can obtain 
through series rearrangement a multi-summation expression for the azimuthal Fourier 
coefficients. Vilenkin's polyspherical coordinates are simply subgroup-type coordinate 
systems which parametrize points on S'^"'^ (for a detailed discussion of these coordinate 
systems see chapter 4 in Cohl (2010) [6]). In this section we will give an explicit example 
of just such an addition theorem on H^. 

The azimuthal Fourier coefficients on H|. expressed in standard hyperspherical 
coordinates ([T]) are given by the functions : [0, oo)^ x [0,7r]^ — )■ R which is defined 
by fl27|) . By expressing P9|) in ci = 3 we obtain 

?/|(x,x') = , ~' 

AirRy sinh r sinh r' 

oo 

X ^(-1)'(2/ + l)P,75'/'+'^ (cosh r<)g;/^+' (cosh r>)P; (cos 7), (59) 

1=0 

1 /2 

where Pi : [—1,1] R is the Legendre polynomial defined by Pi{x) = (a;), or 
through f lSS]) with /i = and v G Nq. By using the addition theorem for hyperspherical 
harmonics (143|) with d = ?> using (cos ^, sin ^ cos 0, sin 6' sin 0) to parametrize points on 
S^, we have since the normalized spherical harmonics are 

the addition theorem for spherical harmonics, namely 

PKCOS7) = |[^^Pr(cos^)Pr(cos^')e™(^-'^'\ (60) 

where cos 7 is given by f l25|) . By combining (1591) and (1601) , reversing the order of the two 
summation symbols, and comparing the result with (!26l) we obtain the following single 
summation addition theorem for the azimuthal Fourier coefficients of a fundamental 
solution of Laplace's equation on H|j, namely since f)^ = AirRH^, 



— ■ °° (I — 

vsmhrsmhr ~~| [L -\- m)\ 



l=\m\ 



X PI'' (cos e)Pr (cos e') ^i72 ^ (cosh r < ) Q^^p (cosh r > ) . 

This addition theorem reduces to the corresponding result ((2.4) in Cohl et al. (2001) 
[9]) in the Euclidean R^ limit as r, r' — )■ O"*". 
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5. Discussion 

Re- arrangement of the multi-summation expressions in section H] is possible through 
modification of the order in which the countably infinite space of quantum numbers is 
summed over in a standard hyperspherical coordinate system, namely 

oo oo I ki fcti-5 'Sd-S 

1=0 K Z=0fci=0fc2=0 fcd_4=0 fcd-3=0 A:d_2=-fcd-3 

oo oo oo oo oo oo 

= E E E ••■EEE- 

fcd_2 = -00 fc£j_3 = |fcd_2l fcd-4=fcd-2 k2=k^ ki=k2 ko=ki 

Similar multi-summation re-arrangements have been accomplished previously for 
azimuthal Fourier coefficients of fundamental solutions for the Laplacian in Euclidean 
space (see for instance Cohl et al. (2000) [11]; Cohl et al. (2001) [9]). Comparison 
of the azimuthal Fourier expansions in section [3] (and in particular ( |36l) ) with re- 
arranged Gegenbauer expansions in section H] (and in particular 0491) ) will yield 
new addition theorems for the special functions representing the azimuthal Fourier 
coefficients of a fundamental solution of the Laplacian on the hyperboloid. These implied 
addition theorems will provide new special function identities for the azimuthal Fourier 
coefficients, which are hyperbolic generalizations of particular associated Legendre 
functions of the second kind. In odd-dimensions, these special functions reduce to 
toroidal harmonics. 
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